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ABSTRACT 


Logistics Operations Command and Control Capability Concept (LOCCC), 
developed by Jeff Grelson in 2000, introduees a new distribution principle to eombat 
elements. This eoneept employs a supporting logisties unit in a general support role and 
eontrols it by a unique eommand eenter in order to minimize the footprint left by 
logisties, improve logistie and taetieal responsiveness, and reduee the “iron mountain” on 
the battlefield. This thesis revisits the mathematieal models and algorithms developed by 
Major Thomas Lenhardt to model LOCCC. We preproeess the network topology in order 
to eonvert it into an equivalent, simplified network that is eomputationally tractable with 
the existing optimization model by using exaet and heuristie algorithms. We show that 
the simplifications and enhancements we propose help us to obtain much faster and better 
quality solutions than using the original, non-simplified networks. For example, in a ten- 
minute run, we ean obtain a solution that is 98% better in some eases. We also apply the 
model to a Turkish Infantry Brigade to evaluate LOCCC with sustainment requirements 
and transportation assets of the Turkish Army. 
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DISCLAIMER 


The reader is warned that all the military terminology used in the Turkish Army 
Manuals and the names of these manuals have been translated by the author for the 
purpose of this thesis. Therefore, it might not be proper to use them in offieial documents. 
All the unit organizations and formations, personnel and equipment quantities, and 
geographical names in the Turkish Army scenario are notional. They have been assumed 
and approximated for this study. Real data is kept classified. The computer code used in 
the thesis cannot be considered validated for use in a different setting than presented in 
this thesis. Any additional application of the source code or the mathematical model is at 
the user’s own risk. 
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EXECUTIVE SUMMARY 


Logistics Operations Command and Control Capability Concept (LOCCC), 
developed by Jeff Grelson in 2000, introduces a new distribution principle to combat 
elements. According to this concept, a Force Service Support Group (FSSG) forms only 
two Mobile Combat Service Support Detachments (MCSSDs) instead of assigning one to 
each Regimental Task Force. One of these MCSSDs supports the entire Ground Combat 
Element in a general support (GS) role while the other adds flexibility to the distribution 
in a direct support role. Jeff Grelson suggests that LOCCC will: (a) help generate a daily 
tasking order after receiving supply requests from supported units at a central command 
and control facility; (b) send the requested support; (c) make necessary modifications to 
the plan as the tactical situation changes; (d) meet immediate supply requests by taking 
advantage of having ready GS units as reserve; and (e) reduce the logistics footprint by 
nearly 60% by lessening the size of the “iron mountain.” 

Major Thomas Lenhardt studied this concept in 2001 and evaluated the existing 
and proposed concepts on how to use the Combat Service Support Element resources of a 
ESSG by solving a vehicle routing problem with demands to be met in specific time 
windows. This thesis revisits the data structures embedded in Major Lenhardfs 
mathematical models and algorithms. By making some simplifications and enhancements 
in these structures, we demonstrate the resulting networks are computationally tractable 
with the same models, while still maintaining all the necessary information required to 
devise an optimal distribution schedule. 

Specific accomplishments of this thesis are: 

• Preprocessing of the network data, which are converted into equivalent data sets that 
are computationally tractable with the existing model formulation. This is 
accomplished by applying two essential steps to the original data: 

o Creating a simplified network: Supply has to be delivered to a small 
number of nodes in the network (the nodes accommodating troops) 
whereas the majority of nodes serve as transition nodes. These transition 
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nodes increase the burden of the model in terms of number of variables 
and equations. We compute all-to-all shortest path distances and then 
extract transition nodes to simplify the network and reduce the model size. 
This simplification does not cause any lack of realism in the model, 
because distances and routes are maintained within our shortest path 
algorithm. Moreover, rounding the travel times is more accurate in the 
simplified network (where we round the total travel time between two 
demand nodes) than in the original network (where we need to round 
times on each leg). 

o Symmetry breaking: The all-to-all connectivity between the demand nodes 
causes redundant arcs. The network is simplified further by extracting 
these replicated arcs, helping the model eliminate redundant solutions. 

These enhancements allow us to obtain better solutions in a reasonable time. For 
example, in one of the cases proposed by Major Lenhardt, the original network would 
provide a solution consisting of 143.85 stons of unmet demand in 10 minutes. After 
our enhancements, we attain a solution of 3.14 stons of unmet demand in the same 
amount of time. The preprocessing time to create the simplified network is, of course, 
negligible. 

• Other managerial and organizational changes in the code, which allow easy access 
and changes to problem data. 

• Application of the model to a Turkish Scenario. 

• Allowing reloading (i.e., redelivery) in the model, in order to enhance its realism. 
This allows us to improve our solution by 2% in some cases of the U.S. Marine 
Corps scenario and to be able to meet the entire demand in the Turkish Army 
scenario. 

• Trying different vehicle combinations in the heuristic algorithm, in order to assess 
how the model solution is affected by these modifications. 



I. INTRODUCTION 


A. BACKGROUND 

A fighting force is composed of two essential elements: war-fighters and 
logisticians. The operational success requires their firm solidarity and harmony, because 
they cannot succeed without each other. History has witnessed many victories and 
breakdowns linked to the success or failure in acquisition and distribution of logistic 
supplies, as well as the heroism or cowardice in conducting tactical operations. 
Napoleon’s and Hitler’s collapse against Russia in different eras, and the United States' 
difficulties to deliver required supplies timely to rapidly advancing forces during 
Operation Iraqi Freedom are examples of how difficult it can be to execute logistics 
efficiently. Success stories also exist in different eras. Conquering commanders have 
always appreciated the role of logistics in combat and have made adequate preparation 
before embarking on a fight. For example, the Ottoman Empire conducted very long 
marches to various areas between the 13* and 18* centuries. Local officials had been 
stocking the necessary supplies and preparing the logistics supply lines on the way to the 
objective months before the army began moving. Similar groundwork was completed by 
General Eisenhower before the Normandy landings and by General Schwarzkopf before 
the Operation Desert Storm. It would be impossible for U.S. troops to succeed in their 
wars across the ocean without proper logistics. Among multiple applications, logistics 
planners can employ optimization models to schedule an optimal distribution of supplies 
to the forces on the battlefield. 

1. Logistic Distribution in a Marine Air Ground Task Force 

A United States Marine Corps (USMC) Marine Air Ground Task Eorce 
(MAGTE) Eorce Service Support Group (ESSG) uses a conventional distribution model, 
which is roughly a refined version of the one used in the World War II era (Grelson 
2000). It consists of a General Support (GS) Combat Service Support Detachment 
(CSSD) and a Direct Support CSSD that assigns its Mobile Combat Service Support 
Detachments (MCSSD) to the supported units (Eigure 1). That model does not require a 
real command and control capability. MCSSD assets essentially receive demands from 
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the units to which they are attached and respond accordingly. MCSSD and supported unit 
commanders have the initiative to decide the vehicle routing schedules. Unfortunately, 
when immediately available reserve assets are limited, headquarters lacks the flexibility 
to take control in case of an urgent demand. Furthermore, transportation assets controlled 
by subordinate MCSSDs are likely to be inefficient because they consider only meeting 
the demands of its attached unit. An emerging concept that keeps transportation assets in 
a general support role is called Logistics Operations Command and Control Capability 
(LOCCC). LOCCC routes vehicles by a unique command center in CSSD to schedule 
deliveries making use of the available loading capacity. Centralized control and 
coordination of the assets also allows the CSSD commander to take immediate actions 
for suddenly rising demands more easily. 



Figure 1. Distribution model used by a MAGTF FSSG (After Gannon, 2000) 


The United States Marine Corps has ten supply classes. Below are short 
definitions of supply classes as defined in Marine Corps Warfighting Publication 
(MCWP) 4-11. 

I. Subsistence including rations, health and welfare items. 

II. Minor end items. 
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III. Petroleum, oil, and lubricants (POLs). 

IV. Construction materials. 

V. Ammunition. 

VI. Personal demand items and nonmilitary sales items. 

VII. Major end items. 

VIII. Medical materiel. 

IX. Repair parts. 

X. Nonmilitary materiel. 

From the above list, water, food, fuel and ammunition expectably occupy most of 
the loading capacity. Accordingly, this thesis takes into account the first, third and fifth 
classes of supplies. Other classes are either negligible in terms of demand or require 
special carrying methods and vehicles; as such, they are not subject to the routing model 
discussed in this thesis. 

Next, we present brief definitions of the service support organizational elements: 

a. Force Service Support Group 

FSSG is the Combat Service Support Element (CSSE) for a Marine 
Expeditionary Eorce (MEE). As defined in MCWP 4-11, it is a grouping of functional 
battalions that provide tactical level ground logistics support to all elements of the MEE. 

b. Combat Service Support Detachment 

CSSDs are special task organizations from various sources tailored to 
meet logistic needs of the units within the structural body of the MAGTE. Its primary 
tasks are to rearm, refuel, and provide limited maintenance repair and supply for a 
supported unit (MCWP 4-24). CSSDs enable the Combat Service Support Commander to 
provide flexible logistic support to different units. CSSDs can be assigned two different 
standard missions. 

(1) Direct Support: Department of Defense Dictionary of Military 
and Associated Terms (JP 1-02) defines direct support as “a mission requiring a force to 
support another specific force and authorizing it to answer directly the supported force’s 
request for assistance.” Supported and supporting units are in a close one-to-one 
relationship. A supported unit requests its requirements directly from supporting units. 
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Higher commanders and headquarters of both sides rarely interfere with the demand and 
supply process. 

(2) General Support: JP 1-02 defines general support as the most 
centralized mission that is given to the supported force as a whole and not to any 
particular subdivision thereof. It gives MAGTF and CSSE commanders the opportunity 
to control the workload and schedule of their subordinates fully and determine the 
priorities of the demands received from the supported units. The proposed LOCCC 
concept revolves around GS. 

c. Mobile Combat Service Support Detachment 

MCSSD is the mobile version of a CSSD with the capability to keep pace 
with a supported maneuver unit. It has the same capabilities as the parent CSSD. 

2, Logistic Distribution in a Turkish Infantry Brigade 

The Turkish Army (TUARM) is executing a transition to a new logistics 
management system. Most of the subordinate units have already started operating in the 
new system while the rest are preparing for the transition. The differences between the 
former system and the new one occur mostly at the strategic level where the automation 
of demand and response processes enable an effective combination of “pull” and “push” 
methods. Almost no changes occur regarding the distribution methods of supplies at the 
brigade level, which is the focus of this thesis. 

The push method is active and based on consumption estimates depending on 
operational tempo. Resources are scheduled to be delivered beforehand according to 
planning factors. The pull method is more “reactive” than the push method, because it 
bases the delivery schedule on actual consumption rates of the units on the field rather 
than on estimates. The supplied units make their requests to the higher commands and 
pull the supply. 

Supply classes are slightly different in TUARM. For example, USMC has ten 
supply classes, whereas TUARM has six, as specified below (TUARM Field Manual 
KKT 100-10, 2003): 
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I. Subsistence and health: Food, fodder (for available animals such as mine 
dogs, security dogs or mules in mountainous regions), water, cigarettes, sanitary 
equipment, medicine, etc. 

II. Main equipment shown in the formation of the unit and spare parts to 
replenish this equipment when necessary. 

III. Petroleum, oil and lubricant (POL). 

IV. Equipment not shown in the formation, but required for special purposes. 

V. Ammunition and explosives. 

VI. Miscellaneous supplies not included in the first five classes, such as 
rescued or captured equipment from enemy hands, psychological warfare equipment, 
maps, water in bulk containers, etc. 

Primarily, TUARM uses four distribution methods to deliver supplies (TUARM 
Field Manual KKT 54-1 1994). 

a. Distribution Methods Used by the Turkish Army 

(1) Take-To-Unit Method: Supplies are delivered to subordinate 
units by the vehicles provided by superior echelons. Generally, class I supplies are 
distributed by this method. In today’s combat operations, the take-to-unit method is 
favored and logistics planners are encouraged to seek ways to exploit this method as 
much as possible. This is because the take-to-unit method allows the higher commanders 
to have a centralized control over transportation assets. The model discussed in this thesis 
purely uses this method. 

(2) Distribution Node Method: Supplies are received by the 
subordinate units from superior units’ depots. Transportation is executed by the vehicles 
of the supported units. Generally class II, IV, and V supplies are distributed by this 
method. 

(3) Mixed Method: This method uses a mixture of the first two 
distribution methods. It is generally used for distribution of class III supplies. 

(4) Vehicle-To-Vehicle Method: Supplies are passed from the 
vehicles of a supporting unit to those of a supported unit at a predetermined meeting 
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point. This method is the least preferred beeause of eongestion and exposure at the 
loading zone. 

b. Distribution Methods for I, III and V Supply Classes 

(1) Class I: Food and water have the majority of weight in this 
supply elass. Supply is delivered by the take-to-unit method down to a brigade eehelon. 
Companies and battalions piek up their shares from the distribution node and take them to 
their kitehens or depots. Sinee this thesis foeuses on brigade level and below, we ean 
aeeommodate this elass of supply in the distribution node method (Figure 2). 



Figure 2. Distribution of class I supplies at brigade level in TUARM 

(2) Class III; Distribution proeess of elass III is very similar to that 
of elass I. Petroleum Oil and Lubrieants (POL) are transported as far forward as possible 
and supporting units use the take-to-unit method as mueh as possible. At the brigade level 
and below, distribution is performed using the distribution node method. Subordinate 
units eollect appropriate empty eontainers and tankers, refill them at the brigade eentral 
distribution node, and deliver them to neeessary vehieles (Figure 3). 
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(3) Class V: The distribution node method is used to deliver 
ammunition to the units. The distribution node is loeated at either army or eorps level. A 
more eentralized distribution model takes plaee for this elass of supply (Figure 4). 
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Figure 4. Distribution of class V supplies at brigade level in TUARM 


In general, TUARM employs the distribution node model for delivery of every 
elass of supply at brigade level and below. Battalions and related eompanies use their 
own vehicles to load supplies from the distribution node of the brigade or the corps. 
Subordinate commanders have the control of the vehicles and they decide which class of 
supply has the priority to use the transportation assets. 


3, Logistics Operations Command and Control Capability Concept 

The LOCCC concept was proposed by Colonel Jeff Grelson (USMC Ret) 
(Grelson 2000). Stated logistic goals of MAGTF are to minimize the footprint left by 
logistics, improve logistic and tactical responsiveness, and reduce the “iron mountain” on 
the battlefield. Those goals require not only high-technology equipment and weapons, but 
also efficient and effective command and control capability. 

Today’s concept of logistics is just a refined version of the one being used in the 
World War II era (Grelson 2000). CSSE forms one CSSD for the general support 

missions of the force, and one MCSSD for each maneuver element. Therefore, the 
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number of logistics units is always one bigger than the number of maneuver units. 
MCSSDs act as if they were attached to their supported units. Such logistics does not 
need a sophisticated command and control capability. Command and control is limited to 
receiving and responding to requests from subordinate units. 

In LOCCC, the CSSE will form only two MCSSDs, “one that will provide 
support to the entire Ground Combat Element in a GS role, and one that will add depth 
and flexibility to the tactical logistics effort in a direct support role” (Grelson 2000). That 
direct support MCSSD does a traveling salesman route (Ahuja et ah, 1993) from the GS 
CSSD to the supported task forces (Gannon 2000) (Eigure 5). Col Grelson suggests that 
EOCCC will: (a) help generate a daily tasking order after receiving supply requests from 
supported units at a central command and control facility; (b) send the requested support; 
(c) make necessary modifications to the plan as the tactical situation changes; (d) meet 
immediate supply requests by taking advantage of having ready GS units as reserve; and 
(e) reduce the logistics footprint by nearly 60% by lessening the size of the “iron 
mountain.” As the size of the iron mountain is reduced, the need for a more effective 
command and control capability increases, because more effort is required to meet 
demands with less supply available. 



Eigure 5. Distribution model proposed in LOCCC concept 
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B, RECENT STUDIES AND THESIS GOAL 

This thesis puts forward an effort to improve the efficiency of the optimization 
model and methods proposed by Major Lenhardt (Lenhardt 2001) to implement the 
LOCCC concept. 

Unlike many other studies on distribution models, Major Lenhardt did not take 
airlift logistics capability into account. Theoretically, aircraft are to be assigned to tactical 
and logistic missions. However, in real life, airlift capability is always exploited by the 
war fighters to be used in tactical jobs. Major Lenhardt included no airlift or sealift 
capability in the model, because he realistically considered that the logistics distribution 
would be assumed by ground transportation assets. His thesis, essentially, evaluated the 
existing and proposed concepts on how to use the CSSE resources of a FSSG by solving 
a vehicle routing problem (VRP) with demands to be met in specific time windows. 
Major Lenhardt coded a discrete event simulation model using the output of an 
optimization model as the input. He concluded that the vehicular speed was the most 
significant component to meet the demands in the required time windows. In order to be 
able to meet the whole demand in the scenario used to evaluate the model, additional 
ground transportation assets and/or aerial logistic support was required. One of the 
downsides of the optimization model he implemented was the computational time 
required to run the full model (i.e., to solve the model exactly, without using heuristics). 
This thesis enhances Major Lenhardt's work to make the models and algorithms used 
computationally more affordable. 

A key enhancement is the way we deal with the network data structure. We 
suggest to convert the given data (e.g., the data originally used by Major Lenhardt) into 
an equivalent data set that is computationally tractable with the original model 
formulation. Other managerial and organizational changes in the code allow easy access 
and changes to problem data. 

This thesis also applies the vehicle routing model used for a USMC MAGTF to a 
Turkish Infantry Brigade. It aims to evaluate the benefits of keeping all the transportation 
assets of the quartermaster company and the ordnance company gathered and directed by 
a central command and control facility. From this facility, supplies are delivered to 


10 



maneuver units on the field. The total number of vehieles in quartermaster and ordnanee 
eompanies in a brigade adds up to 75 during wartime. That number is quite larger than 22 
vehicles available for a CSSD within the body of MAGTF in USMC (Lenhardt 2001). 
This thesis will also explore the extent to which the problem with more vehicles can still 
be solved in affordable time. 

We finally incorporate extra features to the model, such as the ability to reload the 
trucks during the planning time, which is, in fact, a real aspect of the CSSE logistics. 
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II. MODEL FORMULATION 


This chapter presents the mathematical formulation of the VRP whose foundation 
was established by Major Lenhardt (Lenhardt 2001). We discuss: (a) adjustments to 
allow changes in dispatehing order of the vehieles in the heuristic method, (b) 
adjustments in the data to speed up computations and enhance the accuracy of the results, 
and (e) adjustments in the model to incorporate multiple loading of the vehicles. 

We first define and diseuss the VRP specifications, and then describe a 
mathematical model for this problem. 


A. PROBLEM SPECIFICATIONS 

The VRP setting is a distribution network, whieh eomprises nodes (demand sites, 
or points where two or more roads concur) and roads (hereafter ares). Ares are undireeted 
in the network allowing two-way traffic flow, and travel times depend on the arc and 
vehicle type only (eongestion effects, uncertainty due to weather conditions or other 
contingencies are disregarded). 

Nodes have demand only if they accommodate a supported unit. Demand is 
speeified by node, demand type and amount, and a time window when the delivery is 
supposed to be made (otherwise, the non-delivered fraetion of demand is considered 
unmet). In other words, demand is eonsidered met when one or several vehicles distribute 
the required cargo within the specified time window. Our goal is, preeisely, to minimize 
total (weighed) unmet demand. 

Each available vehicle starts and finishes at the origin node, which represents the 
CSSE or the loading depot. Meanwhile, it might be en route from one node to another, 
waiting at a certain node or making delivery. 

Different types of vehicles are employed in the VRP. No maintenance needs or 

attrition are taken into eonsideration. Vehieles have fuel limitations, whieh in turn, limit 

their operating hours and distance per day. All the vehieles in the model maintain a 

constant speed, which is assumed a datum. Delivery times do not depend on the time of 

the day that delivery is made, cargo type or the amount the vehiele is carrying, but solely 
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on vehicle type. Vehicles have maximum capacities for total load and maximum loading 
capacities for each commodity. Commodities are I, III and V USMC classes of supply, 
namely meals ready to eat (MRE), water (bottled and bulk), fuel, and ammunition. 
Certain restrictions apply to transportation of those supply classes. For instance, fuel and 
ammunition cannot be carried together for security reasons. 


B, MATHEMATICAL FORMULATION 

The mathematical model used in this thesis was established by Major Tom A. 
Lenhardt (Lenhardt 2001). The model is restated here with minor changes. 

1. Sets and indices 

T, time periods, tG r = {1,2,...,I r 1} 

C, commodities, ce C 

V, vehicles, vg F 

M, truck types, me M = [HMMWV,FTON,L FS} 

N, nodes, i,je N 

Note: We assume Og A, and node 0 is the origin (CSSE) node. 

A, set of arcs in the network, a = {i,j)e A cz NxN 

2 . Parameters (Data) 

Note: The unit of measure used for commodities is “short ton” (ston) (2,000 
pounds) for the USMC scenario and metric ton (1,000 kilograms) for the TUARM 
scenario. 

typsvm, parameter that takes the value 1 if vehicle v is a truck of type m and 0 
otherwise. 

Note: Each vehicle falls in one and only one type, i.e., E (fPe™ =1, VvG F 

msM 
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deniic^ demand of commodity c at node i (stons) 

travij^ travel time between node i and j through arc a=(ij) for {ij)^ A (number 
of time periods) 

maxTy, maximum route time allowed for vehicle v (number of time periods) 
qvc, capacity of commodity c in vehicle v (stons) 

Note: This parameter is calculated as = Y^type^^q”'\ where q'"‘^ is a given 

m€.M 

capacity of commodity c for vehicles of type m & M 

max q ^, maximum capacity of vehicle v (stons) 

earfyi, earliest delivery time for node i (time period) 

latCi, latest delivery time for node i (time period) 

b.^, unloading time at node i for vehicle v (number of time periods) 

Note: This parameter is calculated as ’ where h™ is a given 

meM 

unloading time at node i for vehicles of type me M 

, penalty per unit for unmet demand of commodity c at node i (penalty 
unit/short ton) 

bigM, big scalar used in calculations for loading the Logistics Vehicle System 
(LVS) and Five Ton (FTON) trucks 

Note: We may set bigM = EE dem^^. 

ieN ceC 

£, small value used in the objective function to discourage vehicles from 
making unnecessary trips. In all examples £ is chosen as 0.00001, which suffices to 
accomplish that goal. 


15 



3, Decision Variables 
Binary Decision Variables 


Xvijt, 1 if vehicle v starts trip through arc A in time period t; 0 otherwise 

Wviu 1 if vehicle v is waiting at i in time period t; 0 otherwise 

Dvit, 1 if vehicle v starts delivering cargo at i in time period t; 0 otherwise 
LWv, 1 if vehicle v acts as a vehicle that transports water 

LFv, 1 if vehicle v acts as a vehicle that transports fuel 

Non-negative Decision Variables 

Svict, quantity of commodity c served by vehicle v at node i in time period t 

(stons) 

Lyc, quantity of commodity c loaded in vehicle v (stons) 

Uic, unmet demand of commodity c at node i (stons) 

4. Mathematical Formulation of the Vehicle Routing Problem 

Minimize 

Z +ZZ Z (1) 

ieN ceC veV {iJ)&A teT\ ieN vg.V t&T\ 

t>latej earlyi<t<latej 


Subject to: 


Z Y,s„„=k.. 

ieN\ teT\ 

dem^^ >0 earlyj <t<late^ 

VvG F,cg C 

(2) 

< dem^^ , Vv g F, / g V, c g C 

dem^^ > 0 , early ^ <t< late^ 

(3) 

Z +Uiy=dem,^, 

teT\ veV 

early 1 <t<latei 

V/ G V, c G C dem^^ > 0 

(4) 

Z^vc ^maxq^, 

ceC 

VvG V 

(5) 
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JeNI 

- ^vi,t+\+ + X^vy,(+1 ’ 

ye VI 

OJ'mA 


VvG E,zg N,t&T\t<\T\ 

(6) 

LW^ + LF^<\, VvgEI 

type..FTON",v = 1 ov type.,Fys”,y = 1 

(V) 

L . +L .. ,.<bisM*LW, 

V, water v, ammo ’ v rt 


Vvg F| 

type.,FTON’\v = 1 or type«Lvs'\v = 1 

(8) 

K"juer^bigM*LF^, 



Vvg F| 

type.,FTON",v = 1 or type„^ys'\v = 1 

(9) 

jeNI 

(0,J)eA 

Vvg F,tG r t > max J[, 

(10) 

iJeNI t^T 

Vvg F 

(11) 

X X <1 

^ vijt ~ ^ vy7,?+;rav^--l “ ’ 

Vvg F, i,j& N,t&T\{i,j),{j,i)^ 

A(n) 

W =1 

Vvg F 

(13) 

o 

II 

Vvg F 

(14) 

K,i = o, 

Vvg F,zg A-{0} 

(15) 


Vvg F,(z,7)g yl 

(16) 

II 

o 

\/v&V,i& N,t&T\t < early, ovt- 

= 1(17) 

Lyc ^ ^VC . 

VvG F,cg C 

(18) 


Vvg V,(l,j)e A,ie N,te T 

(19) 

^vict’^vc’^ic - 0 

Vvg F,zg A,cg C,tG T 

(20) 


5. Modifications to the model for the Turkish Army scenario 

Equations (7) - (9) are slightly revised for the TUARM seenario. The 
eompatibility eonstraints apply to ah vehiele types in this seenario. The modified 
equations are as follows: 

LW^ + LF^< \ , VveV (TUARM-7) 


17 



+ Krammo'' ^ bigM * LW,, VvE V (TUARM-8) 

b^.Vuer ^ bigM *LF^, VveF (TUARM-9) 

6, Description of the Formulation 

a. Objective Function (1) 

We minimize the sum of penalized unmet demand ineurred at all nodes. 
Very small penalties are applied to the vehicles wandering around without delivering 
anything or making notional deliveries. This would be possible otherwise, because there 
might be multiple optimal solutions (e.g., in theory, vehicles could move freely once they 
have served their demand). These penalty values are so small that they do not affect the 
purpose of the objective function while they discourage the vehicles to take needless 
actions. 

b. Loading Constraint (2) 

This constraint prevents the vehicles from delivering more than their 
loads. It also implies that the total quantity of a particular commodity loaded on all of the 
vehicles cannot exceed the total demand for that commodity. 

c. Delivery Constraint (3) 

This constraint prevents the vehicles from delivering more than the node 
demand for every commodity. It also ensures that the delivery is made in the required 
time window, but only when the status of the vehicle is suitable for a delivery (i.e., when 
the vehicle is physically at the incumbent node). 

d. Demand Constraint (4) 

This constraint keeps track of the met and unmet demand for each 
commodity at each node. 

e. Capacity Constraint (5) 

This constraint guarantees that the total load on a vehicle does not exceed 
its maximum capacity. 

f. Balance Constraint (6) 

This constraint ensures that a vehicle is at exactly one state at a particular 
time period. It might be waiting or making delivery at a certain node, or traveling 
between nodes. Trips between nodes are not interrupted. If a vehicle is directed from 
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node to node “ 7 ”, the eompletion of that trip is guaranteed. Figure 6 depiets the 
possible transitions between states at one time step. 



Figure 6 . Description of the vehicle statuses and possible transitions among them 

g. Compatibility Constraints (7) - (9) 

These eonstraints arrange the eompatibility restrietions eoneeming the 
transportation of different types of eommodities on the speeified vehiele types. Fuel 
eannot be earried with any other eommodity exeept MRE. Other eombinations of 
eommodities are allowed to be earried on the vehieles. For the TUARM seenario we 
replaee (7)-(9) by (TUARM-7)-(TUARM-9). 

h. Driver Shift Constraint (10) 

A driver is not allowed to drive more than ten hours a day in both USMC 
and Turkish Army seenarios. A driver is foreed to be at the origin node before he exeeeds 
his maximum driving time. In this formulation, no driver ehange is allowed in the model. 

L Fuel Constraint (11) 

This eonstraint ensures that vehieles do not operate longer than their 
eapabilities based on fuel eapaeities. No refueling and/or reloading is allowed (this 
assumption will be relaxed later in the thesis). Note: Maximum operating hours for all 
types of vehieles in the USMC seenario were ealeulated by Major Lenhardt in his thesis 
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through TMll 240-15/4B. The same ealeulation is obtained from Turkish Army 
Regulation KKY 54-5 for the Turkish Army scenario. 

j. Backtracking Constraint (12) 

This constraint aims to help the solver eliminate unrealistic solutions by 
restricting the vehicles from backtracking on the arcs. 

k. Initial and Final Conditions (13) - (17) 

These constraints ensure all vehicles start and finish at the CSSE node. 

I Maximum Commodity Capacity (18) 

Each commodity can be loaded on a vehicle up to the vehicle’s capacity. 
This constraint keeps the total quantity of a specific commodity on a vehicle below the 
maximum capacity applying to that commodity. 

m. Domains for the Decision Variables (19) -( 20) 

These constraints define the appropriate domains for the binary and non¬ 
negative decision variables. 
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III. SCENARIO DESCRIPTION 


This thesis evaluates two different seenarios. They include the distribution of 
supplies in a Marine Expeditionary Brigade (MEB) of EISMC and in a Mechanized 
Infantry Brigade of TEIARM. The EISMC scenario is directly adopted from Major 
Eenhardt’s thesis and no major changes have been made in order to maintain a fair 
comparison basis after the improvements made in the algorithms, which are described in 
the next chapter. We also attempt to evaluate EOCCC concept in a Turkish territory, with 
Turkish soldiers and transportation assets. 


A. U.S. MARINE CORPS SCENARIO AND DATA 

Three sub-scenarios are run in this scenario group. Essential differences between 
these sub-scenarios are the quantities, time windows and locations of demands. 

The USMC scenario is taken from an exercise at Marine Corps Air Ground 
Combat Center in Twenty-nine Palms, California. Notional data from the Operation Steel 
Knight, the transportation and distribution part of Operation Desert Knight, was used for 
this exercise. 

Sub-scenarios are named as A, B, and C. Scenario A takes account of the full re¬ 
supply of a MEB size unit. Scenario B maintains the same unit locations and time 
windows, but the demands at the nodes are 50% less than the ones in Scenario A in order 
to represent a partial re-supply of the same unit (Eigure 7). 
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Units are indicated in white boxes. The corresponding delivery time windows to 
be met are indicated in gray. Oval tags are the abbreviations for the geographical name of 
the nodes. Square boxes show the time to travel on the associated arc. Time is 
represented in 20-minute time steps, e.g., the 17-21 time window represents a delivery 
requirement between 340 and 420 minutes after the beginning of the operation. This 
network can be used as an overlay to the map of the region, shown in Appendix C. 


Scenario C (Figure 8) represents the distribution of supplies to a smaller special 
task unit. Demands, time windows and the unit locations are completely different from 
the two previous scenarios. 
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Figure 8. Graphical representation of the network for scenario C (After Lenhardt, 

2001 ). 

A Special Task Force is employed in this scenario. Unit types are not given in 
detail in Major Lenhardt’s thesis. 

Sustainment requirements for the units were computed by using a spreadsheet 
model called LOG2000 that was developed by Major Neita Armstrong (Armstrong, 
2000). These figures are deterministic because they depend solely on planning factors. 
Only class I, III, and V supplies are included in the model. Appendix C explains 
computation of sustainment requirements. 

Vehicle types and quantities are from Exercise Desert Knight / Desert Steel 2001, 
Commanding Officer Confirmation Brief (Lenhardt, 2001). There are twenty-two trucks 
divided into three types (Table 1). Maximum operating time for the vehicles is ten hours. 
Delivery times are deterministic and originated by Major Lenhardt from Combat Service 
Support Lield Guide. Vehicles are not subject to limitations on driver availability. 
However, both operating times and fuel capacities are used as constraints in the model. 
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TYPE 

NUMBER 

MAXIMUM 

CAPACITY 

(stons) 

DELIVERY 
TIME (min) 

12.5 TON, TVS 

14 

20.00 

80 

CARGO TRUCK, FTON 

5 

5.00 

40 

UTILITY TRUCK, HMMWV 

3 

1.10 

20 


Table 1. Vehicle characteristics in USMC scenario 


Travel times amid the nodes are proportional to the vehiele speed, whieh is 
assumed eonstant and set to 15 km/h. Travel times on ares are defined in 20-minute time 
steps and rounded up to the nearest integer number (e.g., all of these three time values: 
1.2, 1.7, and 2 would be rounded up to 2 time steps before using them in our model). The 
unit loeations and the time windows in the network are notional as interpreted by Major 
Lenhardt (2001). Some roads that ean never be used in an optimal solution are by-passed 
or ignored to simplify the model. The distribution proeess is demonstrated down to the 
eompany level. In a real eombat operation, platoons and even squads might be loeated far 
from eaeh other and separate deliveries to their loeations might be needed. Companies or 
battalions within the proximity of eaeh other are eonsidered as demand zones. A truek is 
assumed to be at the target even if exaet loeations of subordinate units might be spread 
out. Appendix C shows the summaries of demand zones. 


B, TURKISH ARMY SCENARIO AND DATA 

This seenario is taken from a preparation eourse for Turkish War College. It takes 
plaee at Traee Region in Northwestern Turkey and ineorporates a eorps level unit. We 
have isolated a speeifie segment that eontains the area of responsibility of a meehanized 
infantry brigade. Some adaptations, sueh as appropriately positioning the artillery units 
and defining main supply routes, have been neeessary. Time windows are arranged to 
avoid making the units remain totally off-duty during re-supplying and to make 
subsequent deliveries to neighboring units (Figure 9). The data are notional beeause of 
assumptions and approximations oeeuring. 

Daily sustainment requirements are deterministie and based on planning faetors 
derived from TUARM Direetive KKY 54-5. They eomprise I, III, and V TUARM elass 

supplies. Appendix C explains daily sustainment requirements. 

24 



Vehicle data are all from the war-time records of quartermaster and ordnance 
companies. Some vehicles might have similar capabilities (even if they are of different 
types), so they are given generic names in this scenario. The restrictions on maximum 
operating hours and commodity compatibilities in the USMC scenario apply to this 
scenario, too. Table 2 shows vehicle characteristics. 


TRUCK TYPE 

NUMBER 

MAXIMUM 

CAPACITY 

(tons) 

DELIVERY 
TIME (min) 

LANDROVER 

15 

1.40 

20 

MERCEDES 

37 

2.75 

40 

MAN 

5 

5.00 

60 

TRAILER 

11 

1.50 

20 

TENTON 

2 

10.00 

100 

FUELTANKER SMALL 

3 

10.00 

60 

FUELTANKER BIG 

2 

20.00 

60 

Table 2. Vehicle cl 

liaracteristics in TUARM Scenario 


Travel times are proportional to the speeds of the vehicles. The distribution of 
supplies is preferred to be conducted at night to take advantage of darkness and relax the 
traffic on the roads for operational uses during daytime. In the dark and without 
headlights, trucks are assumed to travel at 16 km/h (KKY 54-5, 1994). Time is 
represented in 20-minute intervals and each travel time is rounded up to the nearest 
integer number. 
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1 Tank Battalion 
Bvigada HQ 
Signal Comp. 


Battalion 
(30-J?) 


1 inf. 
Comp. 

(5-17) 


1 Arty Btr 

(25 - 31) 


1 kif. 
Comp. 

(14-22) 


Comp- 

(24-31) 


Figure 9. Graphical representation of the network for TUARM scenario. 

Oval tags represent all nodes in the network. Units located in certain nodes are 
indicated in white boxes with corresponding delivery time windows in parentheses. 
Square boxes show the time to travel on the associated arc. Time is represented in 20- 
minute time steps. This network can be used as an overlay to the map of the region in 
Appendix C. 
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IV. SOLVING THE VEHICLE ROUTING MODEL 


A. HEURISTIC METHOD 

The VRP model studied in this thesis is hard to solve exaetly, as shown in 
Lenhardt (2001). For example, eaeh USMC scenario has 169,178 equations, 179,185 
continuous variables and 165,412 discrete (binary) variables. Major Lenhardt concluded 
that reaching a guaranteed optimal solution would require an enormous amount of 
computational effort. He used a heuristic technique to reach an acceptable (but sub- 
optimal) solution within affordable time (Table 2). 


scenario 

group size 

computational 

time 

total 

unmet 

demand 

(stons) 

lower 

bound 

(stons) 

linear 

relaxation 

result 

(stons) 

A 

22 (exact method) 

4h 06m 08s 

61.45 

27.56 

7.79 

2 (heuristic) 

39m 18s 

81.86 

N/A 

4 (heuristic) 

Ih 55m 31s 

73.32 

N/A 


22 (exact method) 

4h 05m 18s 

2.87 

0.00 


B 

2 (heuristic) 

13m 32s 

0.44 

N/A 

0.00 


4 (heuristic) 

Hm25s 

0.97 

N/A 



22 (exact method) 

37m 40s 

0.00 

0.00 


C 

2 (heuristic) 

24m 17s 

0.00 

N/A 

0.00 


4 (heuristic) 

2h 09m 00s 

0.00 

N/A 



Table 3. Model results reported in Major Lenhardt’s thesis (Lenhardt, 2001) 


Computational times can be affordable only when a heuristic with small vehicle 
group sizes is used. Heuristics with larger group sizes and the exact method cannot be 
solved within a reasonable time. Major Lenhardt also devised a lower bound value for the 
exact method, which helps us to assess how close we are to the optimal solution. In 
scenario A, the best feasible solution (before the model was interrupted after ten hours, 
where the best solution was found after approximately four hours) was far from the lower 
bound. This yielded an absolute optimality gap of 33.89 stons and a relative optimality 
gap of 55%. For the heuristic, a lower bound can only be obtained by using the linear 
relaxation of the model. For example, in scenario A, we know a solution under 7.79 stons 
of unmet demand can never be achieved. 


Lenhardt's heuristic basically separates the vehicles into groups of a pre¬ 
determined size and approximates the problem solution by using those groups of vehicles 
sequentially. The first group of vehicles distributes the supplies myopically (disregarding 
other vehicles that are available), and turns over the unmet demand to the following 
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group (Figure 10). The proeess repeats until all demand is met or all groups are 
exhausted. 


A trade-off exists between using a heuristie to save eomputational time and 
aehieving an optimal solution. Our heuristie might end up at an optimal or near-optimal 
solution, without being able to prove it. The simplest way to obtain a lower bound for the 
heuristie solution is to solve the linear relaxation of the full model (with all vehieles), but 
in general we expeet this bound to be weak. Configurations of the vehiele groups are 
influential on the heuristie results. We run all the seenarios by using two different 
eonfigurations of the vehieles: groups of two and groups of four vehieles, respeetively. 


Select group size and configuration 


_3 

Select the first of the unr 

r 

outed groups of vehicles 



Solve the problem by considering only the vehicles In the 
current group and dispatch the group 



Leave the unmet demand to the next group of vehicles 


_: 

- Are there any ur 

:_ 

routed groups? 


NO 

Stop and print the results 


Figure 10. Heuristic algorithm (after Lenhardt, 2001) 


We need to introduee some notation in order to deseribe the heuristie algorithm: 
gv, group number for vehiele v 

g, order of the group partieipating in the solution at a heuristie iteration. This 
sealar is used as a eounter, and inereases after eaeh heuristie iteration is finished 
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To formalize the ideas behind the heuristie, let us eonsider the twenty-two 
vehieles in the USMC seenario. Assume that a vehiele group eontains two trueks. Eaeh 
vehiele will have a gv value between 1 and 11 (e.g., gv = 2 for the third and fourth 
vehieles, and gv = 3 for the fifth and sixth). The heuristie algorithm illustrated in Figure 
10 ean be sueeinetly stated in the following steps: 

(1) SetG(e.g., G={l,2,...ll})andgvG G, Vve V. 

(2) Setg:=l. 

(3) Set V’:= { ve F| gv=g}. 

(4) Solve VRP using V’ instead of V. 

(5) If g=\G\, print the solution and STOP. 

(6) Update all demands (subtraeting the demand met by ve F’); inerease g’ by 1 
and return to step 3. 

We also explore modifying the vehiele groups. In his thesis, Major Lenhardt 
ehose to dispateh the vehieles from the first to the last. This grouping ehoiee lets the 
smallest assets deliver their loads earlier and utilizes the large trueks to meet leftover 
demands. We investigate whether reversing the grouping order, whieh allows the larger 
trueks to deliver the bulk of the demand and allows the smaller ones to aoeomplish the 
remaining unmet demand, helps reduee eomputational time or total unmet demand. The 
next ehapter diseusses and eompares the results. 

B, ENHANCEMENTS TO SIMPLIFY THE PROBLEM 

A deeision maker’s priority should be obtaining a reasonable (ideally near- 
optimal) solution within an affordable time. This thesis emphasizes some adjustments and 
improvements to the existing data model in order to make it eomputationally treatable. 

1. Creation of a Simplified Network 

Some nodes in the network provide aeeommodation for military units, so 
requirements and demands for eertain supply elasses arise at those nodes. The remaining 
nodes (with no demand) merely serve for transition purposes. Vehieles either just pass 
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through or wait at those nodes for some time. Transition nodes eonstitute the majority of 
the nodes, thus inereasing the burden of the model in terms of the number of variables 
and equations. We pull those nodes out and simplify the network remarkably by 
following the next steps: 

Step 1: We eompute the shortest paths among all nodes. This is done by using a 
Java (Sun Mierosystems, 2005) eomputer eode implementing the Floyd-Warshall 
algorithm (Abuja et ah, 1993). (See the algorithm in Appendix B.) This algorithm obtains 
a matrix of shortest path distanees and a matrix of predeeessor nodes, whieh enables us to 
maintain the neeessary information about intermediate nodes without keeping them 
explieitly in our VRP model formulation. For example, assume we have A and D demand 
nodes, and B and C transition nodes. Suppose the shortest path between A and D is going 
through B and C. After we remove nodes B and C in the simplified network, the 
predeeessor matrix keeps the travel direetions for getting from A to D (see Figures 11 and 
12 ). 



Figure 11. A sample network before extracting the transition nodes 
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Figure 12. Simplified network 

The network disregards two intermediate nodes and replaees intermediate ares by 
a single are. 


The simplified network has another advantage: Assume in the original network 
transition times on three adjaeent ares are 1.3, 1.3 and 1.3 time periods, respeetively. 
Eaeh of these is rounded as 2, 2 and 2 (beeause our time-phased model requires integral 
travel times), making the total route time equal to 6 time periods (120 minutes), whieh is 
elearly an overestimate of the aetual route time. If we use the aeeurate travel times on 
ares as the input to the Floyd-Warshall algorithm, total travel time for the route would 
add up to 3.9 time periods. This would be rounded to the nearest integer, that is, 4 time 
periods (80 minutes), whieh is a mueh better approximation and gives more flexibility to 
the vehieles to distribute other demands. 

Remark: For eomparison purposes with Major Lenhardt’s results, we still assume 
the rounded values for eaeh are that he used (instead of more aeeurate route travel times). 

Finally, sinee the transition nodes are taken away from the network, vehieles do 
not have an opportunity to wait at those nodes without making any deliveries. In real life, 
this feature will faeilitate the eontrol and eoordination of the vehieles. They will be able 
to wait only at the origin or the demand nodes. 

Step 2: We ereate the simplified network by eliminating all non-demand nodes 
and eonneeting all demand nodes to eaeh other by a direeted are whose length is given by 
the Floyd - Warshall algorithm. 


31 




Consequently, the network in Figure 13 is obtained for scenarios A and B and can 
be compared to the one in Figure 7. 



Figure 13. Network representation for scenarios A and B after node elimination, 

2, Symmetry Breaking 

In spite of a notable reduction in the total number of nodes, the number of arcs is 
still large because of the all-to-all connectivity (Figure 13). Many of those arcs are 
actually passing through the transition nodes that are out of sight here. Therefore, the 
shortest path between two demand nodes (represented as a direct arc) might be of the 
same distance as several arcs in the network. That causes symmetry in the model, i.e., 
multiple solutions in the formulation represent, in practice, the same solution. Symmetry 
can be avoided by extracting the replicated arcs, i.e., those arcs that represent paths 
(between demand nodes) that do not contain intermediate nodes only. This issue can be 
clarified further with the following example; Assume A, D and H are the demand nodes 
in the network. Remaining nodes are transition nodes. The reduced network consists of 
nodes A, D and H only, with arcs (A,D), (D,H) and (A,H). However, if the shortest path 
between A and H uses D, it will not be necessary to consider the arc between A and H 
(the dashed line in Figure 14, whose length would be precisely the sum of distances 
(A,D) plus (D,H)). On the other hand, if the shortest path from A to H does not use D, we 
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need to eonsider a direet are from A to H. This idea will help us eliminate unneeessary 
ares and reduee the size of the problem, whieh is espeeially important in large networks. 



Our symmetry breaking algorithm is depieted in Figure 15. The predeeessor tree 
from the Floyd-Warshall algorithm is used to deeide whether or not every are in the 
simplified network is redundant. 


For all ares {ij) in the simplified network (i.e., network with demand nodes 
only and all-to-all ares) do; 

1. Let k = pred(/,/) 

2. If A: = i, do not remove the are(/,/) from the are list and exit 

Else if A: is a node in the simplified network, remove are(/,/) and exit 
Else set j = k and go baek to step 1. 

Eigure 15. Algorithm to eliminate transition nodes and break symmetry in the 

simplified network, 

Eigure 16 depiets the final network layout after symmetry breaking. Realize the 
enormous simplifieation made with respeet to the original network in Eigure 7. 
Aeeordingly, signifieant reduetion in model size has been aehieved (Eigures 17 and 18). 
This yields eonsiderable time saving in our eomputations. The next ehapter discusses the 
results. 


33 








Figure 16. Network representation for scenarios A and B after breaking the symmetry. 



Figure 17. Network size reduction after the specified improvements in the data 

structures. 
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Figure 18. Reduction in the model size after the specified improvements in the data 

structure. 


C. MODEL ENHANCEMENTS 

The original model allows the vehicles to be loaded at most once during a day. 
Every vehicle leaves the origin node, completes its delivery, and returns to the CSSE 
node; it cannot be reloaded and make a second delivery, regardless how much time is left 
in the day. Eor example, in Eigure 19, we depict a vehicle leaving the origin at time 
period 2, and returning at time period 15. This vehicle waits at the origin, probably 
inefficiently, for the remainder of the day. 

Active time for the vehicle Vehk'le wati ng at the O'lgin njde 

En route or making delivery 


Time 

period 


Eigure 19. Example of working and idle time periods for a vehicle after the first run 

This thesis modifies the implementation of the model to allow vehicles to make a 
second delivery. Conditions required for a second trip and assumptions to conduct it are 
as follows: 

• Part of the total demand cannot be met by dispatching all the vehicles 
once. Unmet demand still exists at some of the nodes. 
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• Vacant times of the vehieles after the first delivery must fit the demand 
time windows of the nodes waiting for a seeond delivery. 

• Reloading times of the vehieles are assumed to be equal to their delivery 
times. 

• When a vehiele is baek to the origin node (garage) from its first delivery, 
it ean get refueled and take a fresh driver, if neeessary. This allows us to 
ignore fuel and erew time eonstraints in the model. 

The working prineiple of this enhaneement is to run the model twiee. Onee the 
model is run for the first time, the aetive and reloading time periods for eaeh vehiele are 
fixed to prevent it from being used again in the seeond run. Then, leftover demands are 
passed to the seeond run (Figure 20). The model seeks an optimal way to exploit the idle 
times of the vehieles (Figure 21). The infiuenee of this improvement is more observable 
in the TUARM seenario due to more suitable time window distributions. In the USMC 
seenario, the majority of the nodes have demands between the time steps 10 and 25. 
Therefore, for example, vehieles with an idle time period between the time steps 25 and 
42 after the first round of delivery eannot be utilized where there might still exists some 
unmet demand to be met between the time steps 10 and 25. The time windows in the 
TUARM seenario is more evenly distributed than those in the USMC seenarios. Thus, the 
vehieles have greater ehanee to find appropriate demand nodes to make their seeond 
deliveries. Results are diseussed in the next ehapter. 
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Figure 20. Algorithm for allowing the vehicles to reload and make a second delivery 



Figure 21. Example of working and idle time periods for a vehicle after the second run 

The vehicle goes for a second delivery during its idle time and meets some 
demand. Model efficiency can be increased by using this approach. 


Remark: An alternative way to handle reloading explicitly in the model would be 
to consider load balance equations for the vehicles at each time period. By adding a sub¬ 
index “t” to the amount of load on the vehicle {Lyct), and a new variable to represent 
loading (or reloading) at the origin {Rvct), we can establish a new formulation with the 


following constraints: 

o 

VI 

VvG V,te T 

(Reload-1) 

^vct ^vc,t-\ ^ ^vict ^vct 

VvG F,cg C,tG T 

(Reload-2) 
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(Reload-3) 


37 














The first constraint would ensure that reloading occurs at the origin node only. 
The second constraint is the load balance equation for each vehicle. The third equation 
ensures the vehicle total capacity is not exceeded at any time. 

Although this formulation requires additional variables instead of since 
these are continuous variables, it is likely that the new model can be solved almost as 
efficiently as the one without reloading. We have not explored this new formulation in 
the scope of this thesis. 
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V. RESULTS 


This chapter presents the results obtained from the exaet and the heuristie 
methods. The USMC seenarios are solved before and after the improvements diseussed in 
Chapter IV. Speeifieally, we foeus our analysis in solution quality and eomputational 
time. Subsequently, TUARM seenario is solved and applieability of LOCCC to a Turkish 
Infantry Brigade is diseussed. All the eomputations are exeeuted on an Intel ® Pentium® 
4 CPU, 2 GHz eomputer with 1 Gb of RAM running under Mierosoft Windows 2000 
operating system. The optimization models are eoded in General Algebraie Modeling 
System (GAMS) (Brooke et al. 1998) and solved by CPLEX 9.0 (ILOG 2004). 

A, U.S. MARINE CORPS SCENARIO 

Major Lenhardt built a model to explore Combat Serviee Support eoneepts. A 
brief summary of his eonelusions is as follows: 

• When a full re-supply is required by a MEB, most of the demand can be 
met by CSSE. If 100% support is required for a MEB, more ground 
transportation assets have to be added to the fleet or airlift has to be 
eonsidered. 

• Partial re-supply of a MEB (50% of full re-supply) ean be handled by 
CSSE. 

• CSSE is eapable of meeting the demands of smaller size units sueh as 
speeial task forees. 

• The VRP model eould be improved by enhaneing model realism, 
deereasing total unmet demand in the solution, and simplifying model 
strueture to ensure a timelier solution. 

Time is one of the most erueial elements in deeision making. Deeision makers are 
usually under time pressure to make eritieal deeisions. Tools that are intended to faeilitate 
the deeision making should be utilized in reasonable time frames. How the node 
elimination and the symmetry breaking lessen the model eomplexity is explained in the 
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following sections. We demonstrate the enhancements in solving times and objective 
function values owing to the simplifications in model structure. The outputs show that a 
noteworthy enhancement is carried out after the improvements are applied to the model. 

As Major Lenhardt noticed, it takes considerably long time to use the exact 
method. Even the heuristic with large vehicle groups requires extensive computational 
effort. Implementing the model in real life might be costly for that reason. (See Table 3 
for results reported in Major Lenhardt’s thesis.) As we will show, the solver makes 
quicker progress in search for an optimal solution after we apply the node elimination and 
the symmetry breaking strategies developed in this thesis. 

Given the computer and the solver version used by Major Lenhardt were not as 
speedy and powerful as the ones used for this thesis, we have executed his computer code 
again using the same resources as for this study, in order to carry out a fair comparison 
and have a better understanding of the actual improvements. 

1. Exact Method 

Time limit is set to 100,000 seconds (27 hrs 47 mins) in order to keep the model 
from running forever. No other restrictions are imposed on CPLEX. Ligures 22 through 
25 show the evolution of the solution value before and after the improvements. Data 
tables for these figures can be found in Appendix D. 

Scenario A is different from the others, because we can prove the demand can 
never be met with the available assets. (A valid lower bound on the optimal solution is 
37.73 stons when the solver is stopped after 27 hours and 47 minutes.) The exact method 
is stopped after 10 hours since there was no further improvement in Major Lenhardt’s 
thesis. It took 4 hours and 6 minutes to reach to 61.45 stons with a relative optimality gap 
of 55%. The single triangle in Ligure 22 represents this solution. The same code has been 
run with a better computer and the most recent solver version. Although initially it 
performs slightly worse than reported by Major Lenhardt, eventually it achieves a better 
objective value. The best feasible integer solution is achieved after 19 hours and 52 
minutes. The execution is stopped at 27 hours and 47 minutes. 
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Either before or after the improvements, the objeetive value evolves negligibly 
after it reaehes the 50 stons mark. Considering this threshold, we observe that, after the 
improvements, we find this solution in 42 minutes using the fully simplified network, 
whereas it takes 10 hours and 36 minutes in the original network. Also from a praetieal 
point of view, if we set the maximum solve time to 10 or 15 minutes (which would be 
acceptable for operational purposes), we can see that the solution after improvements 
(approximately 75 stons of unmet demand) greatly outperforms that before the 
improvements (approximately 300 stons). Figure 23 presents a similar picture of the 
speed gained to close the relative optimality gap. Before the improvements, a relative gap 
under 30% is obtained in 8 hours and 17 minutes, and 10 hours and 42 minutes are 
required for a value below 20%. The same figures are achieved in 42 minutes, and 2 
hours and 33 minutes respectively after the improvements. 



Figure 22. Total unmet demand vs, time before and after improvements (scenario A - 

exact method) 
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Figure 23. Relative gap vs. time before and after improvements (scenario A - exact 

method) 

Gap is computed as ((UB - LB) / UB) where UB (upper bound) is the best 
solution found so far and LB (lower bound) is the best lower bound in the optimal 
solution. 


The benefit of the modeling enhancements developed in this thesis is even more 
clear in scenarios B and C. For scenario B, it takes 4 hours and 5 minutes in Major 
Lenhardt’s implementations in 2001 to meet all the demand. Actually all the demand is 
not met in this scenario, but an unmet quantity of a couple of tons is tolerable for our 
purposes. The same code meets all the demand in 2 hours and 57 minutes with today’s 
computational power. After simplifying the network, it takes only 20 minutes (Figure 
24). 
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Figure 24. Total unmet demand vs. time before and after improvements (scenario B - 

exact method) 

Major Lenhardt’s reports 37 minutes for seenario C. The same eode produeed the 
same optimal solution in 2 hours and 42 minutes today. That differenee is beeause of the 
optimality eriterion defined before the model is started. CPLEX might spend enormous 
time to find a negligibly better solution. However, optimality eriteria eannot be ehanged 
while the model is working. After the improvements, the same objeetive funetion value is 
obtained in 16 minutes (Figure 25). 
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Figure 25. Total unmet demand vs, time before and after improvements (scenario C - 

exact method) 


2, Heuristic Methods 

Heuristics cannot guarantee an optimal solution, but they can be used to find a 
feasible solution to a problem, generally quicker than the exact method. Some limitations 
are applied in order to keep the total solving time reasonably small, which is the primary 
rationale of using heuristics. The model is limited to be solved in 15 minutes with each 
vehicle group, and the absolute and relative optimality gap tolerances are set to be 0.002 
stons and 0.5% respectively. The codes before and after the improvements are executed, 
and the results are illustrated in Figures 26 through 31. Developments in the objective 
value show the transition among the vehicle groups. 

It takes 17 minutes and 35 seconds to obtain the final feasible solution without the 
improvements in the heuristic by using two-vehicle groups. The same case takes only 80 
seconds after the improvements with a slightly better objective function value (Figure 
26). Those times are 48 minutes and 45 minutes respectively for groups of four vehicles 
(Figure 27). It is noteworthy that the model gets to a feasible solution during branch and 
bound (B&B) algorithm more quickly after the improvements are applied. The linear 
Programming (LP) relaxation solution for this scenario is 7.79 stons. This is the only way 
to assess a heuristic solution. Total unmet demand can never go below 7.79 stons for this 
problem. It is not known to us how close to this value the true optimal solution might be. 
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Figure 26. Total unmet demand vs, time before and after improvements (scenario A - 

group size 2) 



Figure 27. Total unmet demand vs, time before and after improvements (scenario A - 

group size 4) 

For scenario B and two-vehicle groups, the solution is obtained in 11 minutes and 
30 seconds before the improvements, whereas it takes only 59 seconds after the 
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Figure 29. Total unmet demand vs, time before and after improvements (scenario B - 

group size 4) 


Similar results ean be seen in Figures 30 and 31 for seenario C. 
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Figure 30. Total unmet demand vs, time before and after improvements (scenario C - 

group size 2) 


UNVETDEMPrO 

(stons) 



Figure 3 1. Total unmet demand vs, time before and after improvements (scenario C - 

group Size 4) 


The solutions ean be improved by using different group arrangements. A very 
large number of combinations are available for vehicle grouping. We only use the one 
that entirely reverses the order. The most notable difference between the two grouping 
orders is the total loading capacities in each group (Table 4). 
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Forward 

Backward 

Vehicles in 
the group 

Total loading 
capacity 
(stons) 

Vehicles in 
the group 

Total loading 
capacity 
(stons) 

1-2-3-4 

8.561 

22-21-20-19 

80.000 

5-6-7-8 

20.000 

18-17-16-15 

80.000 

9-10-11-12 

80.000 

14-13-12-11 

80.000 

13-14-15-16 

80.000 

10-9-8-7 

50.000 

17-18-19-20 

80.000 

6-5-4-3 

16.187 

21-22 

40.000 

2-1 

2.374 


Table 4. Difference between two vehicle grouping orders in terms of total loading 

capacities for groups of 4 vehicles. 


No significant improvement is achieved for scenarios B and C, because they 
already have all the demands met, but we can see the difference in the results obtained 
from scenario A (Table 5). An improvement in the objective function (6.2 stons or 13%) 
emerges as the solving time increases 32 seconds (40%) for two-vehicle groups. For four- 
vehicle groups, however, the improvement gained in the objective value is 4.3 stons 
(9%), which does not seem to offset the extra 1,778 seconds (67%) needed to achieve this 
result. 


Group 

Size 

Grouping 

order 

Total 

Unmet 

Demand 

(stons) 

Resource 

Usage 

Time 

(secs) 

2 

Forward 

52.420 

80.854 

2 

Backward 

46.248 

113.307 

4 

Forward 

49.718 

2722.405 

4 

Backward 

45.428 

4500.000 


Table 5. Differences of the results between grouping orders for scenario A, 


As explained in Section IV.C, an enhancement in the model specification is to 
allow suitable vehicles to make a second delivery after reloading and refueling at the 
origin node. We utilize two-vehicle groups for both deliveries in order to keep the 
computation time minimum, and test scenario A. (Notice the entire demand is met during 
the first delivery in the other two scenarios.) We use the heuristic strategy discussed in 
Section IV.C that allows trucks to make a second delivery, if they have enough idle time 
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at the CSSE node. Vehiele groups are dispatehed in the algorithm from the last to the 
first. The eontribution of the seeond delivery is not very notieeable in this seenario, 
beeause vehieles are busy with the first delivery for most of the day, and their limited idle 
times at the origin node do not fit with required delivery time windows of the units. The 
first round of deliveries takes 113.53 seeonds of solving time and the total unmet demand 
is 46.25 stons. Only one vehiele is able to make a seeond delivery of 1 ston, redueing the 
total unmet demand to 45.25 stons (Figure 32). Other vehieles are not able to meet any 
demand within their idle time frames. The eontribution of seeond deliveries is more 
effeetive in TUARM seenario, as we show in the next seetion. 
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Figure 32. Timeline showing idle and working times of vehicles that are capable of 

making a second delivery in USMC scenario A, 

B, TURKISH ARMY SCENARIO 

This seenario employs more transportation assets than the USMC seenario. That 
means a big swell in the total number of variables and equations. Henee, the model 
beeomes more eomplex, with 379,016 equations and 768,965 deeision variables (even 
after the improvements are applied). Despite the node elimination and the symmetry 
breaking improvements, it still takes an exeeptionally long time to solve the model with 
the exaet method. The model was interrupted at 27 hours and 47 minutes (100,000 
seeonds), but the quantity of delivery was still less than 25% of the total requested 
demand, and the solver was still keeping a zero lower bound. Heuristies of various group 
sizes have been exeeuted. The most myopie heuristie (single vehiele in eaeh group) has 
been tried in addition to two and four-vehiele groups. No eonsiderable differenees exist 
among the objeetive funetion values of different heuristies. However, solution times vary 
substantially (Table 6). The solver is subjeet to the same limitations as in the USMC 
seenario, sueh as a time limit of 15 minutes for eaeh vehiele group, and relative and 
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absolute optimality gaps of 0.5% and 0.02 tons, respectively, for each problem solved in 
the heuristic. 


INITIAL DEMAND = 208.24 stons 

Group 

size 

Grouping 

order 

Computation 
time (secs) 

Unmet 

demand 

(stons) 

1 

Forward 

129 

19.09 

1 

Backward 

110 

18.36 

2 

Forward 

1101 

22.54 

2 

Backward 

822 

17.41 

4 

Forward 

351 

19.32 

4 

Backward 

378 

24.21 


Table 6. Differences of the results among group sizes and orders 


The heuristic with one-vehicle groups provides a satisfactory solution in 
acceptable time; Total unmet demand is less than 9% of total. However, the LP relaxation 
of the full problem is zero, therefore we cannot tell how far these values are from the true 
optimal. Solution times with larger groups increases notably while they produce either 
negligible or no improvements in the total unmet demand. We conclude that a Turkish 
Infantry Brigade can be basically supplied with its own transportation assets (available 
trucks, tankers and trailers in the quartermaster and ordnance companies). The demand 
that is still unmet after the first round of delivery can still be met by allowing reloading of 
some trucks: One-vehicle groups are used in backwards dispatching order. The first 
delivery round is solved in 110 seconds, and the amount of unmet demand turned over to 
the second round of trips is 18.36 tons. We take advantage of idle times as depicted in 
Figure 33. Fourteen vehicles are used for a second delivery, and eventually all the 
demand in the network is met. The reloading capability enables us to meet all the demand 
without using any additional assets. 
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Figure 33. Timeline showing idle and working times of vehicles that are capable of 
reloading and making a second round of trips in TUARM scenario. 
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VI. CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE 

STUDIES 

We adopt Major Lenhardt’s optimization model (Lenhardt, 2001) to represent the 
LOCCC eoneept. We revisit the data structures embedded in his mathematical models 
and algorithms. By making some simplifications and enhancements in these structures, 
we demonstrate the resulting networks are computationally tractable with the existing 
models, while still maintaining all the necessary information required to devise an 
optimal distribution schedule. 

We show that the simplifications and enhancements achieved in this thesis lead to 
a significant improvement in solving time. For example, in a ten-minute run, we can 
obtain a solution that is 98% better in some cases. 

We also test the model to a Turkish Infantry Brigade by using its own sustainment 
requirements and transportation assets. The results show us that the LOCCC concept can 
be applied in this scenario, although the computational effort is substantially higher than 
that of the USMC scenarios. 

Furthermore, we incorporate reloading in the problem, in order to enhance its 
realism. This allows us to improve our solution by 2% in some cases of the USMC 
scenario and meet the entire demand in the Turkish Army scenario. 

Further analyses may improve the model to achieve more realism and reduce its 
complexity. Some possible research opportunities to follow this study are as follows: 

• For a logistics planner, having a delivery at the required time and place is 
essential. In the current model, a time index (representing a 20-minute 
interval) is associated with most variables and constraints, which increases 
the size and difficulty of the model. A possible enhancement would 
consist of monitoring how many times a vehicle delivers before going 
back to the CSSE node. By doing so, time indices can be replaced by 
delivery sequence indices on the variables. Variable numbers would be 
reduced and the model would become less complex. 
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• The CPLEX solver cannot be interrupted to modify the searching criteria 
of the B&B algorithm. For example, optimality criteria have to be initially 
set. B&B sometimes spends an excessive amount of time at a feasible 
solution without improving it. It might be reasonable to accept a “close 
enough” solution to predefined optimality criteria if we have already spent 
much time in improving such solution unsuccessfully. Xpress-Mosel, an 
advanced optimization language developed by Dash Optimization (Dash, 
2005), can be used to overcome this problem by using callbacks that give 
full control over the B&B process. 

• Transportation assets are assumed to have no mechanical failures. Real 
data from previous combats or exercises can be analyzed to get expected 
failure rates for the vehicles. The model can be more realistic if vehicle 
malfunctions are represented. 

• In the current model, the branching strategy assigns higher priorities to 
units with higher demands. A future study can assign priorities according 
to the tactical importance of the units. 

• Time and network resolution can be augmented, e.g., time can be divided 
into smaller intervals of 5 or 10 minutes. 

• The solution of USMC scenario A can be improved and, hopefully, 
finalized. When the solver is interrupted after 100,000 seconds, the lower 
bound is 37.73 stons and the best solution is 45.96 stons. The heuristic 
with four-vehicle groups can already achieve a better solution (45.24 
stons) when the vehicle groups are dispatched backwards. However, the 
true optimal solution to this problem is, as of now, unknown. Similarly, 
we do not know if there is a way to meet the entire demand in the 
TUARM scenario without reloading the vehicles. 

• Reloading and redelivering can be explicitly implemented in the model by 
adding proposed equations in Section IV.C. This would embed the concept 


54 



in the formulation, instead of myopieally solving the model without 
reloading and then allow reloading for idle vehieles. 

• In those eases where entire demand ean be met, there might be multiple 
ways to aehieve it. New goals ean be ineorporated in the model to solve it 
more effeetively. For example, meeting all demand by using the least 
possible number of trueks. 


55 



THIS PAGE INTENTIONALLY LEET BLANK 


56 



APPENDIX A. ACRONYMS AND ABBREVIATIONS 


B&B Branch and Bound 

CSSD Combat Service Support Detaehment 

CSSE Combat Serviee Support Element 

ESSG Eoree Serviee Support Group 

ETON Eive Ton (A vehiele type used in the USMC seenario) 

GAMS General Algebraie Modeling System 

GS General Support 

HMMWV High Mobility Multipurpose Wheeled Vehiele (A vehiele type used in 
the USMC seenario) 

EOCCC Eogistics Operations Command and Control Capability 

EP Einear Programming 

EVS Eogisties Vehiele System (A vehiele type used in the USMC seenario) 

MAGTE Marine Air Ground Task Eoree 

MCSSD Mobile Combat Serviee Support Detaehment 

MCWP Marine Corps Warfighting Publieation 

MEB Marine Expeditionary Brigade 

MEE Marine Expeditionary Eoree 

MRE Meal Ready to Eat 

POE Petroleum, Oil, and Eubrieant 

ston Short Ton 

TUARM Turkish Army 

USMC United States Marine Corps 

VRP Vehiele Routing Problem 
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APPENDIX B. FLOYD - WARSHALL ALGORITHM 


The Floyd-Warshall algorithm (Ahuja et ah, 1993) is a shortest path algorithm to 
find the shortest routes and distanees between every possible node pairs in a network. 
The algorithm produees distanee and predeeessor matriees. The distanee matrix gives us 
the distanee between eaeh node pair. The predeeessor matrix is used to traek the nodes 
along the shortest route between two nodes (from the sink to the souree, traeking 
baekwards). The algorithm below is from Ahuja et al. (1993): 

N: Set of nodes 

A: Set of ares 

pred[i,j]: The node preeeding “j” on the way from “i” to “j” 
d[i,j]: Shortest distanee between nodes “i” and “j” after label eorreetion 
eiji Distanee between nodes “i” and “j” 
n: Cardinality of N 

begin 

for all node pairs [i,j] g NxN do 

d[i,j] : = ooandpred[i,j] := 0; 
for all nodes i e N do d[i,j] := 0; 
for eaeh are (i,j) e A do d[i,j] := eij and pred[i,j] := i; 
for eaeh k := 1 to n do 
for eaeh [i,j] e NxN do 

if d[i,j] > d[i,k] + d[k,j] then 
begin 

d[i,j] :=d[i,k] + d[k,j]; 
pred [i,j] := pred [k,j]; 

end; 

end; 
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APPENDIX C. SCENARIO MAPS AND DATA 


A, U.S. MARINE CORPS SCENARIO 

This scenario is taken from Lenhardt (2001). More elaborate information about 
the data and seenario is available in that study. 

Chapter Ill (Figures 7 and 8) provides the road network with distanees, time 
windows, and unit locations for all three sub-scenarios. Those figures can be used as an 
overlay to the map shown in Figure 34. 
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Figure 34. Map for the USMC scenario (From Lenhardt, 2001) 


Sustainment requirements are derived from LOG2000, whieh is a spreadsheet 
model developed by Neita Armstrong (Armstrong, 2000). Details about the ealeulation of 
sustainment requirements are available from Lenhardt (2001). Demand summaries in 
eaeh seenario are displayed in Tables 7-9. 
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Snpply Class Reqnirements (stons) 

Vicinity of 
Demand Zone 

Units 

# 

PAX 

MRE 

Water 

Fnel 

Ammo 

Sunshine Peak 
(SP3) 

Inf Co 

3 

546 

0.00 

0.00 

0.23 

0.61 



Total: 

546 

0.00 

0.00 

0.23 

0.61 

Sunshine Peak 
(SP4) 

Inf 

HQAVpns 

2 

852 

0.00 

0.00 

5.75 

1.66 



Total: 

852 

0.00 

0.00 

5.75 

1.66 

Lavic Lake (LL4) 

Inf Co 

3 

546 

0.00 

0.00 

0.23 

0.61 



Total: 

546 

0.00 

0.00 

0.23 

0.61 

Gays Pass (GPl) 

Arty Btry 

2 

294 

0.00 

4.75 

23.60 

11.56 



Total: 

294 

0.00 

4.75 

23.60 

11.56 


AAV Co 

3 

588 

1.64 

9.50 

52.18 

1.66 

Emerson Lake 
(ELI) 

AAV HQ 
Co 

1 

368 

1.03 

5.94 

16.71 

5.30 



Total: 

956 

2.67 

15.44 

68.89 

6.96 

Noble Pass (NP2) 

Inf 

HQAVpns 

1 

426 

0.00 

0.00 

2.87 

0.83 



Total: 

426 

0.00 

0.00 

2.87 

0.83 

Black Top (BT2) 

Inf Co 

3 

546 

0.00 

0.00 

0.23 

0.61 



Total: 

546 

0.00 

0.00 

0.23 

0.61 

Noble Pass (NP5) 

Arty Btry 

1 

147 

0.00 

2.37 

11.80 

5.78 



Total: 

147 

0.00 

2.37 

11.80 

5.78 


Tank Co 

3 

258 

0.72 

4.17 

47.48 

1.52 

Black Top (BTl) 

Tank HQ 
Co 

1 

346 

0.96 

5.58 

11.42 

0.56 


LARCo 

1 

139 

0.39 

2.24 

4.26 

1.22 



Total: 

743 

2.07 

11.99 

63.16 

3.30 


Engr Co 

1 

114 

0.32 

1.84 

0.76 

2.28 


Engr HQ 

1 

281 

0.79 

4.54 

25.29 

0.70 

Delta (D4) 

Arty Bn 
HQ 

1 

139 

0.00 

2.25 

6.59 

0.10 


Inf Reg 
HQ 

1 

214 

0.60 

3.46 

6.79 

0.56 



Total: 

748 

1.71 

12.09 

39.43 

3.64 



Grand 

Total: 

5804 

6.45 

46.64 

216.19 

35.56 





Total Demand: 304.84 


Table 7. Summary of sustainment requirements for scenario A (From Lenhardt, 

2001 ). 
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Supply Class Requirements (stuns) 


Vicinity of 
Demand Zone 

Units 

# 

PAX 

MRE 

Water 

Fuel 

Ammo 

Sunshine Peak 
(SP3) 

Inf Co 

3 

546 

0.00 

0.00 

0.12 

0.31 



Total: 

546 

0.00 

0.00 

0.12 

0.31 

Sunshine Peak 
(SP4) 

Inf 

HQAVpns 

2 

852 

0.00 

0.00 

2.88 

0.83 



Total: 

852 

0.00 

0.00 

2.88 

0.83 

Lavic Lake (LL4) 

Inf Co 

3 

546 

0.00 

0.00 

0.12 

0.31 



Total: 

546 

0.00 

0.00 

0.12 

0.31 

Gays Pass (GPl) 

Arty Btry 

2 

294 

0.00 

2.38 

11.80 

5.78 



Total: 

294 

0.00 

2.38 

11.80 

5.78 

Emerson Lake 
(ELI) 

AAV Co 

3 

588 

0.82 

4.75 

26.09 

0.83 

AAV HQ 
Co 

I 

368 

0.52 

2.97 

8.36 

2.65 



Total: 

956 

1.34 

7.72 

34.45 

3.48 

Noble Pass (NP2) 

Inf 

HQAVpns 

I 

426 

0.00 

0.00 

1.44 

0.42 



Total: 

426 

0.00 

0.00 

1.44 

0.42 

Black Top (BT2) 

Inf Co 

3 

546 

0.00 

0.00 

0.12 

0.31 



Total: 

546 

0.00 

0.00 

0.12 

0.31 

Noble Pass (NP5) 

Arty Btry 

I 

147 

0.00 

1.19 

5.90 

2.89 



Total: 

147 

0.00 

1.19 

5.90 

2.89 

Black Top (BTl) 

Tank Co 

3 

258 

0.36 

2.09 

23.74 

0.76 

Tank HQ 
Co 

I 

346 

0.48 

2.79 

5.71 

0.28 

LARCo 

I 

139 

0.20 

1.12 

2.13 

0.61 



Total: 

743 

1.04 

6.00 

31.58 

1.65 

Delta (D4) 

Engr Co 

I 

114 

0.16 

0.92 

0.38 

1.14 

Engr HQ 

I 

281 

0.40 

2.27 

12.65 

0.35 

Arty Bn 
HQ 

I 

139 

0.00 

1.13 

3.30 

0.05 

Inf Reg HQ 

I 

214 

0.30 

1.73 

3.40 

0.28 


Total: 

748 

0.86 

6.05 

19.72 

1.82 

Grand 

Total: 

5804 

3.23 

23.32 

108.095 

17.78 



Total Demand: 152.42 


Table 8. Summary of sustainment requirements for scenario B (From Lenhardt, 

2001 ). 
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Tab 



Supply Class Requirements (stuns) 

Vicinity of 
Demand Zone 

MRE 

Water 

Fuel 

Ammo 

Black Top (BTl) 

1.00 

5.50 

25.00 

2.60 

Delta (Dl) 

0.30 

1.70 

8.00 

0.70 

Noble Pass (NP5) 

1.20 

6.90 

27.00 

3.20 

Lavic Lake (LL4) 

0.40 

2.20 

10.00 

1.00 

Gays Pass (GPl) 

0.50 

2.60 

13.00 

1.20 

Quackenbush Lake 
(QL3) 

0.40 

2.10 

9.00 

0.90 

Emerson Lake (ELI) 

0.50 

2.70 

13.00 

1.25 


Total Demand: 143.85 


e 9. Summary of sustainment requirements for scenario C (From Lenhardt, 


2001 ). 


B. TURKISH ARMY SCENARIO 

The road network for this scenario is depicted in Chapter III. The map below 
shows the actual terrain on which the scenario takes place. It covers Trace Region in 
Northwestern Turkey (Figure 35). 



V-- 

Figure 35. Map for the TUARM scenario 
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The artillery units do not exist in the original seenario. The author adds them to 
taotieally appropriate loeations. The doeuments defining the unit eharaeteristies sueh as 
personnel numbers and tasks, equipment quantities and speeifieations, and unit 
formations are elassified. For this reason, estimates and assumptions are used for those 
eharaeteristies. Logisties Faetors and tables to make the eomputations below are derived 
from the KKT 54-5 Turkish Army Logisties Faetors Manual. 

Food is eonsumed as MREs and weighs 1.91 kg/person/day. Weather eonditions 
are assumed to be normal, so water eonsumption is 53 It/person/day. A liter of water 
weighs 0.982 kg, so total water eonsumed is eomputed as 52.046 kg/person/day. Water 
ean be earried in 5-gallon barrels or tankers tailored for that purpose. Terrain features and 
type of the operation affeet the POL eonsumption rate. The daily POL requirement is 
9.022 It/person/day. Factors of 1.2 for the terrain (moderate smoothness) and 1.5 for 
operation type (defensive) are considered in the calculation. A liter of diesel POL weighs 
0.84 kg, so the total heaviness of PLO is 13.641 kg/person/day. PLO can be carried in 
five or 55-gallon barrels and tankers. Two different planning factors for ammunition 
consumption are available depending on the characteristics of the operation. We assumed 
the average of them for calculation purposes. By doing so, we assume the ammunition 
required is 13.525 kg/person/day. 
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Supply Class Requirements (tons) 

Demand 

Zone 

Units 

# 

PAX 

MRE 

Water 

Fuel 

Ammo 

6 

Arty Batt. 

1 

116 

0.22 

6.04 

1.58 

1.57 

11 

Arty Batt. 

2 

232 

0.44 

12.07 

3.16 

3.14 

14 

Tank Co. 

4 

432 

0.83 

22.48 

5.89 

5.84 

16 

Tank Co. 

4 

432 

0.83 

22.48 

5.89 

5.84 

HQ Co. 

1 

224 

0.43 

11.66 

3.06 

3.03 

Signal Co. 

1 

120 

0.23 

6.25 

1.64 

1.62 

Engineer 

Co. 

1 

183 

0.35 

9.52 

2.50 

2.48 

TOTAL: 

959 

1.83 

1.83 

49.91 

13.08 

26 

Infantry Co. 

1 

153 

0.29 

7.96 

2.09 

2.07 

27 

Tank Co. 

1 

108 

0.21 

5.62 

1.47 

1.46 

28 

Infantry Co. 

1 

153 

0.29 

7.96 

2.09 

2.07 

30 

Infantry Co. 

1 

153 

0.29 

7.96 

2.09 

2.07 

31 

Infantry Co. 

1 

153 

0.29 

7.96 

2.09 

2.07 

32 

Tank Co. 

1 

108 

0.21 

5.62 

1.47 

1.46 


Grand Total: 

2567 

4.90 

4.90 

133.60 

35.02 


Table 10. Summary of sustainment requirements for the TUARM scenario. 
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APPENDIX D. COMPUTATIONAL RESULTS 


Computational results are graphieally displayed in Chapter V. The tables 
presented in this appendix provide the numerieal basis for those graphs. For the exaet 
methods, the only bound imposed on the model is the solving time (set to 100,000 
seeonds). For heuristie implementations, the maximum solution time is 15 minutes at 
every heuristie iteration and the absolute and relative optimality gap toleranees are set to 
be 0.002 stons and 0.5%, respeetively. 


A, U.S. MARINE CORPS SCENARIO 


Before improvements 

After improvements 

Time 

(secs) 

Best integer 
(stons) 

Best node 
(stons) 

Gap 

Time 

(secs) 

Best integer 
(stons) 

Best node 
(stons) 

Gap 

0 

304.84 

11.01 

1.00 

0 

304.84 

11.01 

1.00 

2051 

178.33 

37.73 

0.79 

314 

105.22 

37.73 

0.64 

3162 

177.48 

37.73 

0.79 

336 

98.99 

37.73 

0.62 

3301 

172.71 

37.73 

0.78 

557 

87.04 

37.73 

0.57 

3440 

121.73 

37.73 

0.69 

601 

75.29 

37.73 

0.50 

3579 

101.05 

37.73 

0.63 

623 

74.68 

37.73 

0.49 

3995 

100.98 

37.73 

0.63 

645 

69.48 

37.73 

0.46 

4273 

95.98 

37.73 

0.61 

690 

69.21 

37.73 

0.45 

4412 

85.60 

37.73 

0.56 

734 

68.75 

37.73 

0.45 

4968 

82.21 

37.73 

0.54 

745 

68.73 

37.73 

0.45 

5107 

81.83 

37.73 

0.54 

889 

67.84 

37.73 

0.44 

5523 

81.38 

37.73 

0.54 

1508 

65.71 

37.73 

0.43 

5801 

81.00 

37.73 

0.53 

1795 

62.02 

37.73 

0.39 

6079 

78.64 

37.73 

0.52 

2524 

48.80 

37.73 

0.23 

6218 

77.00 

37.73 

0.51 

4068 

48.80 

37.73 

0.23 

6287 

74.42 

37.73 

0.49 

5880 

48.80 

37.73 

0.23 

14273 

73.73 

37.73 

0.49 

7542 

48.80 

37.73 

0.23 

15523 

73.07 

37.73 

0.48 

8257 

48.62 

37.73 

0.22 

15759 

72.46 

37.73 

0.48 

9179 

*45.95 

37.73 

0.18 

15940 

57.87 

37.73 

0.35 


29829 

51.26 

37.73 

0.26 
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Before improvements 

After improvements 

Time 

(secs) 

Best integer 
(stons) 

Best node 
(stons) 

Gap 

Time Best integer Best node „ 
(secs) (stons) (stons) ^ 

38162 

47.33 

37.73 

0.20 


38579 

46.72 

37.73 

0.19 

38857 

46.64 

37.73 

0.19 

43718 

44.00 

37.73 

0.14 

71495 

*43.96 

37.73 

0.14 


Table 11. Results of scenario A (exact method). 


* These figures are the best feasible solutions before the run is interrupted after 100,000 
seeonds (27 hours and 47 minutes). No further improvement is aehieved. 


Before improvements 

After improvements 

Vehicle 

Cumulative 

Total unmet 
demand 
(stons) 

Vehicle 

Cumulative 

Total unmet 
demand 
(stons) 

group 

time (secs) 

group 

time (secs) 

VI -V4 

35 

297.27 

VI -V4 

4 

296.08 

V5-V8 

78 

277.27 

V5-V8 

14 

276.08 

V9-V12 

978 

210.38 

V9-V12 

914 

206.76 

V13-V16 

1878 

147.38 

V13-V16 

1814 

143.76 

V17-V20 

2778 

84.38 

V17-V20 

2714 

80.76 

V21 -V22 

2874 

52.88 

V21 -V22 

2722 

49.72 


Table 12. Results of scenario A (heuristic with four-vehicle groups). 


Table 13. 


Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

V19-V22 

900 

228.40 

V15 - V18 

1800 

165.40 

Vll - V14 

2700 

102.40 

V7-VIO 

3600 

60.90 

V3 - V7 

4500 

47.80 

VI -V2 

4501 

45.43 


Results of scenario A (heuristic with four-vehic 


e groups) (vehicles 


dispatched backwards). 
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Before improvements 

After improvements 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

VI -V2 

10 

303.45 

VI -V2 

1 

302.27 

V3-V4 

20 

298.59 

V3-V4 

2 

297.27 

V5-V6 

33 

288.59 

V5-V6 

4 

287.27 

V7-V8 

40 

278.59 

V7-V8 

6 

277.27 

V9-V10 

280 

240.77 

V9-V10 

48 

242.05 

Vll -V12 

427 

209.27 

Vll -V12 

55 

209.92 

V13-V14 

538 

177.77 

V13-V14 

64 

178.42 

V15-V16 

602 

146.27 

V15-V16 

67 

146.92 

V17-V18 

728 

114.77 

V17-V18 

72 

115.42 

V19-V20 

890 

83.27 

V19-V20 

77 

83.92 

V21 - V22 

1055 

53,31 

V21 - V22 

81 

52,42 


Table 14. Results of scenario A (heuristic with two-vehicle groups). 


Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

V21 -V22 

6 

264.64 

V19-V20 

11 

231.07 

V17-V18 

33 

199.57 

V15-V16 

47 

168.07 

V13-V14 

69 

136.57 

Vll -V12 

80 

105.07 

V9-V10 

86 

73.57 

V7-V8 

87 

63.57 

V5-V6 

89 

53.57 

V3-V4 

112 

48.62 

VI -V2 

113 

46,25 


Table 15. Results of scenario A (heuristic with two-vehicle groups) (vehicles dispatched 

backwards), 
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B 

iefore improvements 

After improvements 

Time (secs) 

Best integer 
(stons) 

Best node 
(stons) 

Time (secs) 

Best integer 
(stons) 

Best node 
(stons) 

0 

152.42 

0.00 

0 

152.42 

0.00 

1530 

89.84 

0.00 

232 

54.98 

0.00 

1657 

64.69 

0.00 

296 

31.98 

0.00 

1911 

63.92 

0.00 

680 

23.85 

0.00 

3308 

54.75 

0.00 

701 

22.11 

0.00 

3944 

53.16 

0.00 

723 

12.49 

0.00 

4071 

48.51 

0.00 

744 

10.60 

0.00 

4198 

46.57 

0.00 

787 

8.62 

0.00 

4325 

35.08 

0.00 

802 

8.10 

0.00 

5214 

30.90 

0.00 

1171 

4.40 

0.00 

6865 

30.73 

0.00 

1218 

0.24 

0.00 

7246 

22.59 

0.00 

1226 

0.20 

0.00 

7284 

22.47 

0.00 

2366 

0.20 

0.00 

8644 

13.93 

0.00 


9152 

9.74 

0.00 

9660 

8.67 

0.00 

12455 

8.30 

0.00 

12709 

4.04 

0.00 

13064 

3.31 

0.00 

14233 

0.32 

0.00 

17477 

0.20 

0.00 

23024 

0.20 

0.00 

31802 

0.12 

0.00 

37538 

0.12 

0.00 


Table 16. Results of scenario B (exact method). 


Before improvement 

After improvement 

Vehicle 

Cumulative 

Total unmet 
demand 
(stons) 

Vehicle 

Cumulative 

Total unmet 
demand 
(stons) 

group 

time (secs) 

group 

time (secs) 

VI -V4 

30 

143.95 

VI -V4 

3 

143.95 

V5-V8 

930 

124.59 

V5-V8 

903 

125.41 

V9-V12 

1830 

60.55 

V9-V12 

1803 

61.37 

V13-V16 

2730 

14.06 

V13-V16 

2403 

14.88 

V17-V20 

3630 

1.70 

V17-V20 

3303 

1.91 
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Before improvement 

After improvement 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

V21 - V22 

3662 

0.97 

V21 -V22 

3306 

0.77 


Table 17. Results of scenario B (heuristic with four-vehicle groups). 


Vehicle group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

V19-V22 

900 

85.99 

V15 - V18 

1800 

27.54 

Vll - V14 

2047 

5.10 

V7-VIO 

2625 

0.96 

V3 - V7 

3525 

0.63 

VI -V2 

3526 

0.32 


Table 18. Results of scenario B (heuristic with four-vehicle groups) (vehicles 

dispatched backwards). 


Before improvement 

After improvement 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

VI -V2 

11 

150.14 

VI -V2 

1 

150.14 

V3-V4 

25 

143.95 

V3-V4 

3 

143.95 

V5-V6 

40 

133.95 

V5-V6 

5 

133.95 

V7-V8 

66 

124.56 

V7-V8 

8 

124.56 

V9-V10 

122 

92.02 

V9-V10 

17 

92.02 

Vll -V12 

265 

60.52 

Vll -V12 

28 

60.52 

V13-V14 

408 

29.02 

V13-V14 

35 

29.02 

V15-V16 

429 

14.03 

V15-V16 

36 

14.03 

V17-V18 

572 

6.56 

V17-V18 

49 

6.56 

V19-V20 

644 

2.33 

V19-V20 

56 

2.33 

V21 - V22 

691 

1.08 

V21-V22 

59 

1.08 


Table 19. Results of scenario B (heuristic with two-vehicle groups). 
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Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

V21 -V22 

57 

118.53 

V19-V20 

61 

85.99 

V17-V18 

78 

54.49 

V15-V16 

95 

27.54 

V13-V14 

97 

12.55 

Vll -V12 

103 

5.10 

V9-V10 

107 

1.39 

V7-V8 

117 

0.96 

V5-V6 

137 

0.94 

V3-V4 

140 

0.63 

VI -V2 

141 

0.32 


Table 20. Results of scenario B (heuristic with two-vehicle groups) (vehicles dispatched 

backwards). 


B 

iefore improvements 

After improvements 

Time (secs) 

Best integer 
(stons) 

Best node 
(stons) 

Time 

(secs) 

Best integer 
(stons) 

Best node 
(stons) 

0 

143.85 

0.00 

0 

143.85 

0.00 

1213 

56.45 

0.00 

133 

15.90 

0.00 

1808 

49.30 

0.00 

315 

7.00 

0.00 

2328 

45.15 

0.00 

340 

3.20 

0.00 

2477 

38.75 

0.00 

567 

3.14 

0.00 

2849 

34.55 

0.00 

962 

0.00 

0.00 

3518 

29.41 

0.00 


3592 

28.23 

0.00 

3667 

22.18 

0.00 

3741 

20.99 

0.00 

3964 

11.79 

0.00 

4931 

10.00 

0.00 

7681 

8.75 

0.00 

7756 

4.95 

0.00 

8648 

2.40 

0.00 

9786 

0.35 

0.00 

10842 

0.00 

0.00 

Tal 

ble 21. Results of scenario C (exact method). 
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Before improvements 

After improvements 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

> 

1 

39 

135.29 

VI -V4 

5 

135.29 

V5-V8 

939 

116.78 

V5-V8 

905 

116.82 

V9-V12 

1839 

52.78 

V9-V12 

1805 

52.32 

V13-V16 

2739 

13.22 

V13-V16 

2705 

11.97 

V17-V20 

3639 

o 

o 

V17-V20 

2712 

0.01 

V21 -V22 

3642 

0.00 

V21 -V22 

2712 

0.00 


Table 22. Results of scenario C (heuristic with four-vehicle groups). 


Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

V19-V22 

900 

77.45 

V15 - V18 

1800 

30.70 

Vll - V14 

2417 

5.60 

V7-VIO 

2420 

0.00 

V3 - V7 

2421 

0.00 

VI -V2 

2421 

0.00 


Table 23. Results of scenario C (heuristic with four-vehicle groups) (vehicles 

dispatched backwards). 


Before improvements 

After improvement 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

(N 

> 

1 

15 

141.48 

VI -V2 

1 

141.48 

V3-V4 

27 

135.29 

V3-V4 

2 

135.29 

V5-V6 

46 

125.84 

V5-V6 

5 

125.84 

V7-V8 

946 

118.19 

V7-V8 

65 

116.54 
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Before improvements 

After improvement 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

V9-V10 

1046 

85.69 

V9-V10 

82 

84.04 

Vll -V12 

1215 

54.19 

Vll -V12 

no 

52.54 

V13-V14 

1346 

31.19 

V13-V14 

121 

29.14 

V15-V16 

1459 

13.84 

V15-V16 

126 

10.14 

V17-V18 

1499 

2.41 

V17-V18 

129 

1.21 

V19-V20 

1505 

o 

o 

o 

V19-V20 

129 

0.00 

V21 -V22 

1506 

o 

o 

o 

V21 -V22 

129 

0.00 


Table 24. 


Results of scenario C (heuristic with two-vehicle groups). 


Vehicle 

group 

Cumulative 
time (secs) 

Total unmet 
demand 
(stons) 

V21 - V22 

38 

109.75 

V19-V20 

87 

77.45 

V17-V18 

109 

50.85 

V15-V16 

143 

30.70 

V13-V14 

144 

13.35 

Vll -V12 

147 

5.60 

V9-V10 

148 

0.00 

V7-V8 

148 

0.00 

V5-V6 

148 

0.00 

V3-V4 

149 

0.00 

VI -V2 

149 

0.00 


Table 25. Results of scenario C (heuristic with two-vehicle groups) (vehicles dispatched 

backwards). 


B, TURKISH ARMY SCENARIO 

Three different vehiele grouping sizes have been used in the heuristies. The 
results are given in Chapter V (Table 6). 
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